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Population and host plant differentiation of the 
sweetpotato whitefly, Bemisia tabaci ( Homoptera: 
Aleyrodidae ) , in East, South and Southwest China 
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Abstract: The sweetpotato whitefly, Bemisia tabaci ( Gennadius ) , is a species complex that is considered an 
important agricultural pest because of its diversity in biotypes, capacity to transmit viruses, high resistance to 
chemicals and rapid physiological development. Its biotype status in Jiangsu, Guangdong and Yunnan provinces 
was determined by using mitochondrial DNA COI markers and correlated with their host plants and locations. 
Results indicated that, in Guangdong and Yunnan, the local unresolved biotypes and the already known biotypes 
of B. tabaci (B and Q) coexist, the former is relatively more dependent on its host plants than the invasive 
biotypes. No Q biotype of B. tabaci has been found in Guangdong and no local biotype has been found in Jiangsu 
in this study. Multitudinous host plants and the inherent polyphagous characteristics enhance the wide invasion of 
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B and Q biotypes in different counties. 
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1 INTRODUCTION 


whitefly 
( Gennadius) (Hemiptera; Aleyrodidae ) , also known 


The  sweetpotato Bemisia tabaci 
as the cotton or tobacco whitefly, has become one of 


the most devastating tropical and sub-tropical 


worldwide agricultural pests ( Byrne and Bellows, 
1991) in the past 20 years ( Anthony et al., 1995; De 
Barro, 1995; Perring, 2001). This cosmopolitan pest 
affects the yield of a broad range of agricultural, fiber, 
vegetable and ornamental crops (Cahill et al., 1996) 
by causing damage through feeding directly on the 
sap, and through 
honeydews that produces sooty mould. These sooty 
mold fungi, using honeydew as a substrate, colonize 
further 
photosynthesis and ultimately resulting in reduced 
quality of fruit and fiber ( Byrne and Bellows, 1991). 
In addition, B. tabaci is the vector of several 
economically important plant viral-pathogens, these 


plant indirectly excretion of 


contaminated surfaces, interfering with 


could be more than 110 plant viruses, most of them 
belong to begomovirus ( Geminiviridae ), which has 
been recognized as the most important emerging plant 


virus group in sub-tropical and tropical world regions 
(Brown, 1990, 2000; Martin et al., 2000; Jones, 
2003 ; Mughra et al., 2008). B. tabaci is a genetically 
diverse group (Frohlich et al., 1999; De Barro et al., 
2000) and its genetic complexity was first recognized 
in the 1950s when morphologically indistinguishable 
populations were reported to differ in host range, host 
plant adaptability, and plant virus-transmission 
capabilities ( Bird, 1957; Mound, 1963; Costa and 
Russell, 1975; Bird and Maramorosch, 1978). This 
led to the development of the concept that B. tabaci 
was composed of a series of biotypes ( Costa and 
Brown, 1991; Bedford et al., 1994). Since then, 
there has been a proliferation of biotypes described with 
the current count exceeding 20 ( Perring, 2001) and 
more recently, its global populations known have been 
grouped into 12 genetically distinct groups ( Boykin et 
al., 2007). The considerable genetic diversity of B. 
tabaci, while 
characters, with various complexity of their biological 


retaining identical morphological 


interactions and their varying abilities to interbreed, 
have led to the belief that B. 
cryptic species complex ( Bedford et al., 1994; Brown 
et al., 1995; Byrne et al., 1995; Moya et al., 2001; 


tabaci represents a 
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Pascual and Callejas, 2004; De Barro et al., 2006; 
Qiu et al., 2007). 


B. tabaci outbreaks caused by native and 
introduced variants are often not immediately 
recognized until populations have reached 


uncontrollable levels. Although not all populations 
have the same propensity to achieve pest or vector 
status and not all of them are capable of long-distance 
spread, the B and Q biotypes are exceptions known to 
be invasive and have spread worldwide ( Frohlich 
et al., 1999; De Barro et al., 2000, 2005; De Barro, 
2005; Martinez-Carrillo and Brown 2007; Qiu et al., 
2007). Furthermore, in B. tabaci biotypes, B and Q 
are polyphagous, and other biotypes might differ from 
them in host plant specialization. It has already been 
reported that the B. tabaci species complex has over 
600 host plants, including numerous economically 
significant horticultural and ornamental species 
(Mound and Halsey, 1978; Secker et al., 1998). The 
reasons behind, the worldwide spread of B and Q 
biotypes are complicated, but it has been thought that, 
it could be associated with the international trades of 
ornamental and other economically important crops 
( Brown et al., 1995; Perring, 2001; Reitz, 2007). 
In China, B. 


serious pest until its outbreak in the mid-1990s, which 


tabaci was not considered as a 


was due to the invasion of the B biotype ( De Barro et 
al., 2005). This biotype has caused serious damage to 
crops, and 
ornamental plants in more than 30 provinces in China 
(Qiu et al., 2007, 2008). Moreover, the Q biotype, 
which was recorded in 2004 in Yunnan Province, 
Southwest China in 2003 (Chu et al., 2005) , has also 
been found on various host plants in Beijing ( Zhang et 
al., 2005), Jiangsu (Shen, 2007), and Zhejiang 
provinces ( Xu et al., 2006). In addition to the 
invasive biotypes, some native populations or 
haplotypes of B. tabaci, such as the Cv and ZhJ-2 
populations, also cause serious damage to their specific 
host plants (Qiu et al., 2006; Zang and Liu, 2007). 


Preliminary field investigation indicated that most of the 


a wide range of vegetables, fiber 


native B. tabaci populations are oligophagous. 
It has always been problematic to identify these 


of B. 


molecular 


various biotypes tabaci morphologically, 


however, sequence analysis of the 
mitochondrial cytochrome oxidase I gene (mtCOI) has 
been recognized as one of the most useful molecular 
markers for differentiating B. tabaci populations and 
for identifying these genetically diverse populations 
(Bird and Maramorosch, 1978; Costa and Brown, 
1991; Brown et al., 1995; Guirao et al., 1997; 
Frohlich et al., 1999; Berry et al., 2004; Brown, 
2000; Qiu et al., 2007). Population differentiation 


plays an important role in revealing an organism’ s 
evolutionary history, including the B. tabaci species 


complex, but this has not been clearly reported in 
China. In this study, the biotype/haplotype and host 
plants of B. tabaci populations in eastern, southern 
and southwestern parts of China are investigated by 
sequencing the mitochondrial COI gene of various B. 
tabaci populations, collected from multiple plant hosts 
throughout China; the results may show us the 
relationships between population differentiation, host 
plants and geographical distribution of the B. tabaci 
samples, which will help us to estimate the spread of 
different B. tabaci populations in the future. 


2 MATERIALS AND METHODS 


2.1 Insects 

Whitefly adults were collected from the fields 
using a hand-held aspirator and preserved immediately 
in 95% ethanol for later processing, and their host 
plants were taken back to the laboratory for species 
identification simultaneously. A total of 167 samples, 
each approximately 100 - 500 individuals, were 
collected from 11 different provinces during 2006 - 
2008. B. tabaci populations collected from Jiangsu, 
Guangdong and Yunnan provinces were investigated as 
representatives in this study due to the rich 
biodiversities of both B. tabaci populations and their 
host plants, all the B. tabaci samples from these 
provinces are abundant ( 350 - 450 adults) in 
population and their host plants are from various 
species. Their locations and host plants were showed in 
Table 1. 
2.2 Biotype identification of the B. tabaci samples 
2.2.1 DNA extraction: DNA was extracted from 
individual insect of each sample. The B. 
samples were firstly washed using double distilled water 


tabaci 


to remove alcohol, then they were homogenized 
individually in 200 wL lysis buffer ( 1% SDS, 
10 mmol/L Tris-HCl, pH 8.0, 25 mmol/L EDTA, 25 
mmol/L NaCl, 1.0 mg/mL proteinase K) for 30 s, in 
a 1.5 mL microcentrifuge tube. The homogenate was 
incubated at 55°C for 2 -— 3 h, after which the 
homogenate was centrifuged for 5 min at 5 000 rpm and 
the suspension was moved to a new 1. 5 mL 
microcentrifuge tube and used for PCR amplification 
directly or stored at —20°C for further use. 

2.2.2 PCR amplification: The mitochondrial COI 
gene was used to the biotype of the 
representative samples of B. tabaci. Primers used to 
amplify the mtCOI gene fragment were C1-J-2195 (5’ 
TTGATTTTTTGGTCATCCAGAAGT-3’) and L2-N- 
3014 ( 5’-TCCAATGCACTAATCTGCCATATTA-3’). 
PCR reactions (25 uL) contained 2.0 uL of each 
template DNA, 2 pmol/L MgCl, , 25 pmol/L for each 
primer, 2.5 mmol/L of each dNTP, 1.0 uL Express 
Taq polymerase (5 U/ uL, TaKaRa) and 1 x PCR 


identify 
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buffer. PCR amplifications were done in a Bio-Rad 
thermocycler and included a precycle denaturation step 
for 2 min at 94°C ; followed by 35 cycles at 95°C for 1 
min, 50°C for 1 min, and 72°C for 1 min; and a final 
extension at 72°C for 10 min. The PCR products were 


electrophoresised on 0. 8% agarose gels, and bands 
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an ultraviolet light. Those products that yielded a band 
with the expected size of 820 — 850 bp were sequenced 
on an ABI377 automated DNA Sequencer ( PE Applied 
Biosystems). At least 5 individuals of each B. tabaci 
samples were repeated for PCR application and COI 


were visualized using ethidium bromide staining under 


sequencing. 


Table 1 Samples collected from different locations and host plants in China 


Codes 


Collection sites 


Ast-GGD Zhaoqing, Guangdong 
Con-GGD1 Yingde, Guangdong 

Con-GGD2 Zhanjiang, Guangdong 
Con-GGD3 Yangjiang, Guangdong 


Con-GGD4 Shaoguan, Guangdong 


Biotype 


Non-B/Q 
Non-B/Q 
Non-B/Q 
Non-B/Q 
Non-B/Q 


Con-GGD5 Guangzhou, Guangdong Non-B/Q 














Eup-GGD Shenzhen, Guangdong 
Leg-GGD Meizhou, Guangdong 
Sol-YUN1 Baoshan, Yunnan 
Sol-YUN2  Mengla, Yunnan 
Sol-YUN3 Yuanjiang, Yunnan 
Sol-YUN4 Jinghong, Yunnan 
Sol-YUNS Yuanjiang, Yunnan 
Sol-YUN6  lincang, Yunnan 
Con-YUN1 Jinghong, Yunnan 
Con-YUN2 Ruili, Yunnan 
Con-YUN4 Menghai, Yunnan 
Con-YUNS  Pingbian, Yunnan 
Eup-YUN Kunming, Yunnan 
Sol-JSU1 Dongtai, Jiangsu 
Sol-JSU2 Suqian, Jiangsu 
Sol-JSU3 Nantong, Jiangsu 
Sol-JSU4 Haimen, Jiangsu 
Con-JSU1 Changshu, Jiangsu 
Con-JSU2 Dafeng Jiangsu 
Cru-JSU Dafeng, Jiangsu 
Cuc-JSU Dafeng, Jiangsu 
Mal-JSU Xinghua, Jiangsu 
Vit-JSU Dafeng, Jiangsu 


2.3 Data analysis 


For the 5 repeated COI sequences of each B. 
tabact population, 


Non-B/Q 
Non-B/Q 
Non-B/Q 
Non-B/Q 





Non-B/Q 
Q 
Non-B/Q 
Non-B/Q 





Non-B/Q 
B 
Non-B/Q 
Non-B/Q 

Q 


B 
Q 
B 
Q 
Q 
Q 
B 
Q 
B 
Q 


they were firstly analyzed and 
aligned with DNAStar ( LasergeneS. 0), and the most 
accurate sequence was selected to represent this 











GenBank 
Host plant fost Pant accession Locations 
Y no. of mt COI 

Lactuca sativa ( lettuce) Asteraceae  FU192034 23° 02' 53” N, 112° 27/54” E 
Ipomoea batatas (sweet potato) Convolvulaceae EU192036 24° 00’ 09” N, 113° 27’ 20” E 
Ipomoea batatas( sweet potato) Convolvulaceae EU192037 21° 16’ 14” N, 110° 21’ 33” E 
Ipomoea batatas (sweet potato) Convolvulaceae EU192038 21° 51’ 27” N, 111° 58°57" E 
Ipomoea batatas (sweet potato) Convolvulaceae EU192040 24° 48’ 39” N, 113° 35' 60” E 
Ipomoea batatas (sweet potato) Convolvulaceae EU192042 23° 07’ 44"N, 113° 15/51” E 
Codiaeum variegatum (croton) Euphorbiaceae EU192043 22° 32’ 36"N, 114° 03’ 28” E 
Vicia faba (horsebean) Leguminosae EU192039 24° 17’ 18” N, 116° 07’ 20" E 
Lycopersicon esculentum (tomato) Solanaceae EU192052 25° 06’ 43” N, 99° 09’ 41” E 
Lycopersicon esculentum (tomato) Solanaceae EU192054 21° 28’ 00” N, 101° 35’ 00” E 
Solanum melongena (eggplant ) Solanaceae EU192060 23° 48’ 20" N, 102° 05’ 35” E 
Lycopersicon esculentum (tomato) Solanaceae EU192061 22° 00’ 45” N, 100° 47’ 36” E 
Solanum melongena (eggplant ) Solanaceae EU192060 23° 48’ 20" N, 102° 05’ 35” E 
Lycopersicon esculentum (tomato) Solanaceae EU192056 22° 31’ 60” N, 99° 55' 60” E 
Ipomoea batatas (sweet potato) Convolvulaceae EU192053 22° 00’ 45” N, 100° 47’ 36” E 
Ipomoea batatas (sweet potato) Convolvulaceae EU192055 24° 00’ 589” N, 97° 51' 11” E 
Ipomoea batatas (sweet potato) Convolvulaceae EU192058 21° 5858” N, 100° 25’ 43”F 
Ipomoea batatas (sweet potato) Convolvulaceae EU192059 22° 59’ 06” N, 103° 41’ 10” E 
Euphorbia pulcherrima (poinsettia) Euphorbiaceae EU192049 25° 02' 20” N, 102° 43’ 19” E 
Solanum melongena (eggplant ) Solanaceae EU192070 35° 50' 45” N, 120° 18’ 05” E 
Solanum melongena (eggplant) Solanaceae EU192071 33° 56' 57” N, 118° 17’ 45” E 
Solanum melongena (eggplant ) Solanaceae EU263628 32° 10' 14" N, 121° 04’ 28” E 
Lycopersicon esculentum (tomato) Solanaceae EU263629 31° 53’ 32" N, 121° 09 12” E 
Lablab purpureus (hyacinth-bean ) Convolvulaceae EU192072 31° 38’ 39” N, 120° 44' 17" E 
Ipomoea purpurea (sweet potato) Convolvulaceae EU263630 33° 12’ 00” N, 120° 28’ 00" E 
Brassica napus (rape) Cruciferae EU263625 33° 12’ 00” N, 120° 28’ 00” E 
Cucumis sativus (cucumber) Cucurbitaceae EU263633 33° 12’ 00” N, 120° 28’ 00" E 
Gossypium herbaceum ( cotton) Malvaceae EU263627 32° 36' 25" N, 119° 00’ 25” E 
Parthenocissus tricuspidata( Boston ivy) Vitaceae EU263631 33° 12’ 00" N, 120° 28’ 00” E 





population. Then, the representative sequences from 
different B. tabaci populations were compared using 
ClustalX1. 83. The genetic distances between different 
biotypes of B. tabaci were calculated based on their 


COI sequences, and their phylogenetic relationships 


10 #4 


were calculated using ME ( Minimum Evolution ) 
method in MEGA v.3.1 (Bootstrap = 1 000). Some 
COI sequences of B. tabaci populations from Spain, 
Morocco and USA as well as from China were 
downloaded from GenBank and used as reference 
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sequences. The A biotype from USA (DQ174542) and 
the greenhouse whitefly Trialeurodes vaporariorum, 
( AF418672 ) were used as a species and genus 
outgroup in the phylogenetic tree (Table 2). 


Table 2 Samples downloaded from GenBank and used as reference in phylogenetic tree 


Location 


Code/ biotype 


Morocco/Q N/A 


Nador, Morocco 


Spain/Q Spain N/A 


Host plant 


Host plant family GenBank accession no. of mt CO1 


N/A AJ517769 


N/A EU427719 


Eup-USA /B 
Fab-GGD/B 
Con-Zhj-2-1 
Con-Zhj-2-2 
USA/A 


Trialeurodes vaporariorum 


N/A: not available. 


Florida, USA 
Guangzhou, China 
Zhejiang, China 
Zhejiang, China 
USA 


Bangalore, India 


Euphorbia pulcherrima (poinsettia ) 
Glycine max (soybean ) 
lpomoea purpurea ( sweet potato ) 
lpomoea purpurea ( sweet potato ) 
N/A 
Phyllanthus emblica ( Amla) 


66 pm :.Con-YUNS5 


Euphorbiaceae 
Fabaceae 
Convolvulaceae 
Convolvulaceae 
N/A 


Phyllanthaceae 


EU 





EU 


427727 
192032 


DQ309077 


AJ867557 


DQ174542 


AF418672 
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Strict consensus phylogenetic tree of Bemisia tabaci based on CO I sequences 


The tree has been constructed by the ME ( Minimum Evolution) method in MEGA v. 3.1 (Bootstrap = 1 000) to show relationships among different 
Bemisia tabaci biotypes. The sequence of the greenhouse whitefly Trialeurodes vaporariorum ( AF418672 TV) was used as outgroup. 


1136 Be et Acta Entomologica Sinica 


3 RESULTS 


3.1 Phylogenetic relationships of different B. 
tabaci populations 

The phylogenetic relationship of all the B. tabaci 
populations collected from Guangdong, Jiangsu and 
Yunnan provinces is shown in Fig. 1. Results revealed 
that all the Chinese populations are clustered into two 
groups, the non-B/Q group and the B/Q group. For the 
non-B/Q group, there are two branches (branch I and 
IT) ; samples in the first branch are from Yunnan, and 
all of their host plants belonged to the Convolvulaceae 
family, the biotype of this branch was named Yur-1 
(Yunnan unresolved-1 ) in this study. Meanwhile, 
samples in the second branch were divided into two sub- 
branches, the samples in Il-a sub-branch are also from 
Yunnan, but their host plants belong to the Solanaceae 
family, the biotype of this sub-branch was named as 
Yur-2 temporarily in present study. The H-b sub-branch 
includes several populations from Guangdong; their host 
plants are more diverse, including Convolvulaceae, 
Asteraceae, Leguminosae, Euphorbiaceae. The biotype 
of these samples includes the already known biotypes 
Zhj-2 , Cv as well as one unresolved biotype, which was 
named Gur ( Guangdong 


unresolved ) biotype 


52 4 


temporarily in this study, and we also found that the 
Zhj-2 biotype was more close to the Gur population, 
than the Cv biotype. 

Phylogenetic relationships of the B. 
populations in the B and Q group are simpler than that of 
non-B/Q group. All the Q biotype samples from Yunnan, 


tabaci 


Jiangsu and two reference populations from Morocco and 
Spain were clustered into the III branch, and all the B 
biotype samples were clustered into the IV branch. The 
host plants of the Q and B biotype populations are more 
varied, including seven plant families, Convolvulaceae, 
Euphorbiaceae, Solanaceae, Cruciferae, Cucurbitaceae , 
Malvaceae and Vitaceae. 
3.2 Genetic distances between different biotypes 
of B. tabaci Chinese populations 

The genetic distances between different biotypes of 
B. tabaci are shown in Table 3. Results indicated that 
the genetic distances of Cv with Zhj-2 and Gur are 
0.0988 and 0. 966, respectively, while that between 
Zhj-2 and Gur is 0. 0019, which indicates the closer 
relationship of Zhj-2 and Gur. However, the genetic 
distances between different branches are waved due to 
population differentiation. For example, the distance 
between B and Q is very close (0. 0608 ) while that 
between Yur-1 and Yur-2 is much larger (0. 1662). 


Table 3 Genetic distances among different biotypes’ genetic groups of whiteflies based on their CO I sequences 


Biotype/genetic groups B Q 
B Q 0. 0608 
Cv 0. 1832 0. 1793 
Zhj-2 0.1712 0.1770 0. 0988 
Gur 0.1717 0. 1765 0. 0966 
Yur-l 0. 1685 0. 1746 0. 1490 
Yur-2 0. 1832 0. 1789 0. 1098 


* Biotype/haplotype that is unresolved in Guangdong and Yunnan provinces. 


4 DISCUSSION 


The first report of B. tabaci in China was in 1949 
(Chou, 1949), and its biotype diversity has not been 
considered until the invasion of the exotic B biotype 


(Qiu et al., 2005). In China, 
Guangdong and Jiangsu, it has been reported that there 


especially in 


are at least 200 plant species that can be attacked by 
B. tabaci (Qiu et al., 2001; Zhou et al., 2004). On 
the other hand, many local or unresolved biotypes/ 
haplotypes of B. tabaci populations have been found 
recently in different limited locations, and the known 
diversity of this pest is still increasing. Boykin et al. 
(2007 ) resolved the world phylogeny of B. tabaci 
species complex into 12 major genetic groups. The two 


Zhj-2 Gur * Yur-1 * Yur-2 * 
0.0019 
0. 1536 0.1717 
0.1145 0.1144 0. 1662 


most invasive B. tabaci, B and Q biotypes were 
classified into the Mediterranean/Asia Minor/ Africa 
group and Mediterranean group respectively while other 
all Non B/Q biotypes as the China are classified in 
China group by Boykin et al. (2007). 

The B. tabaci samples collected from Guangdong, 
Jiangsu and Yunnan provinces were analyzed in this 
study; it has been found that there are B, Q and at 
least two local unresolved biotypes ( Yur-1 and Yur-2 ) 
in Yunnan. In Guangdong, there are B, Cv, Zhyj-2 
and a local unresolved biotype (Gur) while Q biotype 
found in the 


Guangdong. However, in Jiangsu province, only the 


was not samples collected from 


two invasive biotypes, B and Q were found in this 


study. B. 


tabaci is considered a globally quarantine 
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species in regards to international trade in a range of 
plant and plant product commodities, and it also has 
been suggested to use strict quarantine measurement 
during trade inside the country to stop the movement of 
these biotypes (both native and invasive biotypes). 

In the case of host plant specialization, it is 
interesting that all the non-B/Q biotype populations in 
Yunnan province are divided into two groups; Yur-1 
and Yur-2 and their host plants are quite simple; all of 
Yur-1 belonged to Convolvulaceae while those of Yur-2 
belonged to Solanaceae family. However, those of the 
non-B/Q biotypes from Guangdong province are quite 
complex (Fig. 1). It may indicate that the non-B/Q 
biotype populations in Yunnan province are more 
dependent on their host plants, compared to those B. 
tabaci populations from Guangdong. Geographic 
barriers alone could not explain all the observed genetic 
variations while differential host plant utilization 
explained one of the reasons for genetic variation in 
non-B/Q samples but could not explain the overall 
cause for these variations. Furthermore, roles of mating 
interference and Wolbachia in developing the genetic 
structure could not be ignored ( De Barro, 2005; 
Ahmed, 2008). Chinese populations of B. tabaci used 
in this study, have already been investigated to be 
positive for Wolbachia infections (Ahmed, 2008). As 
far as Q and B biotype B. tabaci concerned, they had 
been collected from various host plants which are 
divided into seven plant families. This might confirm 
the hypothesis that a wide range of host plants is one of 
the reasons that the exotic B and Q biotype populations 
could settle down easily in new places through avoiding 
population from dying out ( Brown et al., 1995 ). 
Furthermore, no Q biotype has been found in 
Guangdong till now, although we do not know if it 
already exists. Therefore, more monitoring work should 
be done in the future. 

In conclusion, the results of this study suggest 
that, in eastern, southern and southwestern China, the 
local unresolved biotypes and the already known 
biotypes of B. tabaci (B and Q) coexist and the non- 
B/Q biotypes relatively more dependent on their host 
plants than the invasive biotypes. More work should 
been done to investigate and control the concurrent 
damages from both invasive and indigenous biotypes of 
B. tabaci in the future. 
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